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Abstract 

■ Neglecting the effect of particle production at the moment of bubble nucleation, the spectrum 

. of created particles during the bubble expansion is evaluated in the thin-wall approximation. It 

CN \ is shown that the expanding thin- walled bubble makes the dominant contribution to the particle 

K»" ' production. 
O " 



In this paper, as a step toward full understanding of the matter, we extend our analysis [1] and 
^ ' investigate the particle creation during the process of bubble growth, i.e. after the bubble nucleation, 
in the thin wall approximation. It has been shown that in the thin-wall approximation the particle 
^ ' production at the moment of bubble nucleation is strongly suppressed [1]. On the other hand, 
^ . numerical calculations show that the particle production by the thick-walled bubble is very intensive 
^ ! [2]. The latter observation has led us to consider the particle creation which occurs during the process 
^ ' of bubble expansion in the thin-wall approximation. The question we ask is the following: can one 
T— I ' neglect in the thin wall case the particle production after the bubble nucleation in comparison with 
2 ■ the number of particles produced during the nucleation process? The expanding bubble is described 
. by the 0(3,1) invariant function which is an analytic continuation of the bounce solution [3]. In 
4^ ', the thin wall case, the bounce looks like a large four- dimensional spherical bubble with a thin wall 
CL( separating the false vacuum without from the true vacuum within. As the bubble expands, the wall 
^ traces out the hyperboloid. The bounce solution provides the classical background with respect to 
which the quantum fluctuations are defined. The second quantization picture for the fluctuation 
'-j . field in process of false vacuum decay through the barrier penetration is thoroughly examined in 
^ ', [4, 5]. Throughout this paper the signature for the Minkowskian space-time metric is (—,+,+,+) 
■ and h = c = 1 is assumed. For the calculation of particle spectrum we use the standard formalism 
developed in [4, 5]. 

The lagrangian density is 

C = -\{d,cl^f-U{ct^l (1) 

where 

f/(0) = ^<f{<p - 2af - e{<f/2a^ - 30Vl6a^), (2) 

and a. A, e are positive parameters. The potential (2) has a local minimum at 0/ = 0, f/(0) = 0, a 
global minimum at 0f = 2a, U{2a) = — e, and a local maximum at 4>top = 8Aa^/ (8Aa^ + 3e), U {(ptop) = 
128A^a^^(2Aa'^ + e)/(8Aa^ + 3e)'^. The validity condition for the thin-wall approximation is a^A/e ^ 1, 
see [3]. Assuming the thin wall approximation the bounce solution with a good accuracy may be 
written as 

0fe = a{l — tanh(/i[p — i?])}, (3) 

where /i = a\/X and the bubble radius R = 2/i^/eA. Following the standard formalism it is convenient 
to use the coordinates which respect the symmetry of the background solution (ph, i.e. 0(4) and 
0(3, 1) in Euclidean and Minkowskian regions respectively. The coordinate system in the Euclidean 



region is given by (p, X, 0, ip) where {9, ip) are usual angle coordinates on two-dimensional sphere and 
(p, x) are related to r = |a;| and r by 



r = psin(x), r = -pcos(x), 
< X < 7r/2, < p < oo. 

The coordinates in the Minkowskian region are obtained by the replacement (p, x) 
{-it -iXu), which yields 

r = ^sinh(xAf), t = ^ cosh(xAf), 
< Xm < oo, < ^ < oo. 



(4) 



The equation governing the fluctuation field $ reads 



$ = 0, 



(5) 



(6) 



where denotes the partial derivative with respect to ^ and L is the Laplacian operator on three- 
dimensional unit hyperboloid. Expanding $ in terms of harmonic functions on the three-dimensional 
hyperboloid Ypi^, —LF'Ypim = (1 +P^)i^zm) the Eq.(6) for the mode function takes the form 



(p2 + 1/4) 



/ = o. 



(7) 



The analytic continuation of / to the Euclidean region is performed by setting ^ ^ ip. Under this 
continuation one obtains the equation for the mode function in the under-barrier region 



5 = 0, 



(8) 



where g is related to / by the asymptotic boundary condition f{^) = g{—iC) at ^ ^ 0. At r — > — oo 
the field is in false vacuum state and correspondingly the fluctuation field satisfies the vanishing 
boundary condition when p — > oo. Since ^/"(^b) 4a^A when p — > oo the asymptotic solution to 
Eq. (8) satisfying the vanishing boundary condition is given by 

g{p) = a{p)^Kip{2aVXp), 

where Kip{x) is the modified Bessel function of the second kind. On the other hand, the asymptotic 
solution to Eq.(7) when ^ — > oo will have the form 



where H^^ {x) and Hl^ [x] are the Hankel functions of the first and second kinds, respectively. The 
spectrum of created particles n{p) is given by 

n{p) -■ 



r{2) 



1 



Cl(p)/C2(p)p - 1| 



(9) 



For detailed description of the general formalism concerning the problem see [4, 5]. In what follows, 
we shall assume that the particle production does not occur at the moment of bubble nucleation. 



To judge the accuracy of this assumption we first evaluate the spectrum of cretaed particles at the 
moment of bubble nucleation. Using the bounce solution (3) for [/"{(ph) one obtains 

U"{(f)b) = 2/x2{3tanh^(//[p - R]) - 1} + e-term. 

For the sake of simplicity we approximate the well of U"{4>b) in the region \p — R\ < by the 
square well 

where for a reasonable approximation we assume 1.5 < d <2. We neglect the contribution coming 
from the e-term because at = this term equals zero and at ~ a this is of order eA///^ ^ //^. 
The solution of Eq. (8) is given by 



9{p) = < 



aiy/plip{2iip) + a2y/pKip {2iip), p < 

as^Hi^\Vdpp) + a.^HffiVdpp), < p < R+, (H) 
^Kip{2pp), p > R+, 



where R± = R± /i^^, overbar denotes complex conjugation, H-p\x) = cxp{p'n-)H-p\x), and Iip{x) 
is the modified Bessel function of the first kind. Since in the thin-wall approximation pR± 3> 1 for 
matching of g at the junction points R± one can use the asymptotic formulas. For x ^ 1 and <^ x 
the asymptotic expressions are given by [6] 

\/2'nx ^ 2.x ^ \ t:x 

For the coefficients ai, 02, 03, 04 one obtains 

aa oc (2 - i\fd) e-'(^^^R+-^/^\ 04 oc -(2 -t- tVd) e^(^/^^+-V4)^ 
where the proportionality coefficients are the same for 03, 04 and 

ai oc V7F(4 + rf)sin(2v^)e"^^^-, 02 oc {(4 - rf) sin(2v^) + 4A/dcos(2v^)} e^^^-/v^, 

the proportionality coefficients are the same for ai, 02- Then for small values of momentum the 
number of created particles n{p) at the moment of bubble nucleation is evaluated as 

For large values of momentum h{p) decreases at least as [1] 

n{p) ~ exp(— 2p7r). (14) 

For evaluating the spectrum of created particles during the bubble expansion one has to solve 
Eq.(7) with the asymptotic boundary condition /(^) = «0 at ^ — > 0. The assumption that the 
particle production docs not occur at the moment of bubble nucleation implies that the solution to 
Eq.(7) in the region ^ < i?_ is given by [4] 



Under this assumption the solution to Eq.(7) reads 



/(O 



Ci^flLip{^/d|^i) + Ci^Kip{\fd^ii), R- < ^ < R+, 



(15) 



where Lip[x) = ^{Iip{x) + I-ip{x)}. Both Kip{x) and Lip{x) are real valued and even functions of p. 
Any of Lip{x), Kip{x) and Iip{x) can be constructed from the remaining two functions by the identity 

■^ipi^) — Lip{x) — i— — ^ — ^-Kip(^x). 

TT 

Using the asymptotic formulas (12) the matching of / at the junction points R± gives 

C3 oc —exp{—^/d|J,R_){^/d + 2i)^/^T, C4 oc exp(\/d/ii?-)(2i — \/d)/\/7r, 

where the proportionality coefficients are the same for C3, C4 and 

ci oc i4\/dcosh(2\/d) + (d - 4) sinh(2\/d), C2 oc -{d + 4) sinh(2\/d), 

again the proportionality coefficients are the same for ci, C2. Thus the spectrum of created particles 
during the bubble expansion for small values of momentum is approximately given by 



n(p) ~ ^^^-^smh.'^(2Vd). 
16a 

For large values of p and x one can use the following asymptotic formulae [6] 



(16) 



Kip{x) ~ 



27re-f^/2 



Lip[x) 



(p2 -x2)V4 



27r(p2-x2)V4 



cos(q;(x) — 7r/4) + — - sin(ci;(x) — 7r/4) I , 
12p I 



sin(Q;(a;) — 7r/4) — — - cos(Q;(a;) — 7r/4) 
12p 



(17) 



gi(/3(a;)-7r/4) I ^ 



12p 



7r(p2+x2)V4 

where Q;(a;) = pln{{p+{p'^—x'^)^/^)/x} — {p'^—x'^y/'^ and /3(a;) = pln{a;/(p+(p2+a;2)^/2)} + (p2_^2;2^i/2_ 
Assume p 3> /xi? then the matching at the junction point R- gives 



C3 OC (27r)V2e-f-/2 



C4 oc -(27r)-V'e^'^/2 



y/djiR. 



^ [sin (7) — icos(7)] — ^ 



p 



-[cos(7) + isin(7)] + 



2p 
HR- 



Vdsm('y) H — j= cos(7) 

yd 

-^sin(7) — \/dcos(7) 



VduR-" '"" ■ ' ■ 2p 

where the proportionality coefficients are the same for C3, C4, and 7 = a{VdijLR-) — 7r/4. Using these 
coefficients and asymptotic formulae (17) from matching at R+ one obtains that |ci/c2| behaves as 
p'^/{liR)^. Correspondingly, for large values of momentum 



(18) 



Thus, in contrast to Eq.(14) the particle production by the expanding thin- walled bubble follows 
a power law in p when p — > oo. From Eq.(18) one concludes that the contribution to the number 
of created particles coming from the ultraviolet region is finite. Prom Eq.(16) one sees that even 
in the thin-wall case the particle creation is not really suppressed if the expansion of the bubble 
is considered. However, in this case the particle production is not so intensive as in the thick- wall 
approximation [2]. This result is not strictly correct since we neglected the constraint on the particle 
production which arises due to fact that the proper fluctuation fleld associated with the tunneling 
is the transverse part of total fluctuation fleld with respect to the bounce solution [7] . Summarizing 
the above results one concludes that in the thin-wall approximation particle production is dominated 
by the expanding bubble. 
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